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Abstract: We extend our earlier results delineating the supersymmetry (SUSY) reach 
of the CERN Large Hadron Collider operating at a centre-of-mass energy yfs = 7 TeV 
to integrated luminosities in the range 5-30 fb" 1 . Our results are presented within the 
paradigm minimal supergravity model (mSUGRA or CMSSM). Using a 6-dimensional grid 
of cuts for the optimization of signal to background ratio - including missing Et~ we find 
for nig ~ nig an LHC 5a SUSY discovery reach of rrig ~ 1.3, 1.4, 1.5 and 1.6 TeV for 
5, 10, 20 and 30 fb -1 , respectively. For rriq 3> m§, the corresponding reach is instead 
m g ~ 0.8, 0.9, 1.0 and 1.05 TeV, for the same integrated luminosities. 
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1. Introduction 

In 2011, the CERN Large Hadron Collider has produced proton-proton collisions at a 
centre-of-mass energy y/s = 7 TeV (LHC7) and has enabled both ATLAS and CMS ex- 
periments to each accumulate over 5 fb _1 of useful data. The current plan is to resume 
running with pp collisions in early 2012, with a goal to amass in the vicinity of 10-30 fb -1 
of usable data. The 2012 run will likely be followed by a shut down for ~ 2.5 years so that 
various upgrades may be implemented; turn-on at or near design energy of y/s = 14 TeV 
is then expected around 2015. 

While many LHC analyses are focused on the elusive Higgs boson, the search for weak 
scale supersymmetry (SUSY) [jD remains an important part of the LHC program. In a 
previous paper ||], we presented projections for the LHC7 5<r discovery reach for SUSY 
in the paradigm minimal supergravity (mSUGRA or CMSSM) model[||]. In that study, 
we presented discovery strategies for early SUSY discovery and made projections for the 
LHC7 reach for a variety of integrated luminosities ranging from 100 pb" 1 up to 2 fb -1 , 
well beyond what was then expected to be delivered in the entire 7 TeV run. LHC reach 
projections for yfs = 14 TeV (LHC14) have been reported in earlier studies 

Recent analyses (performed within the mSUGRA model) of the LHC7 data by the 
ATLAS H and CMS|| experiments based on just ~ 1 fb -1 of integrated luminosity have 
found no indication of SUSY so far, yielding 95% CL lower limits of rriq ~ m g ~ 1 TeV for 
comparable gluino and squark masses, and m g ~ 0.6 TeV for the case where rriq 3> rrig. It 
is worth emphasizing that although all squarks are by assumption degenerate within the 
mSUGRA framework, the squark mass limit cited above arises mostly from signals for first 
generation squarks that are much more copiously produced from qq and qg initial states 
than their second and third generation cousins. In other words, the LHC7 squark limit 
really applies to up and down type squarks - other squark flavors may be significantly 
lighter than the quoted bounds. These LHC7 bounds do not apply to third generation 
squarks or to electroweak-inos, the only sparticles with significant couplings to the Higgs 
sector and to which the naturalness arguments that yield upper mass bounds on sparticles 
apply. Indeed models with 0(10 — 100) TeV gluinos and first generation sfermions but with 
sub-TeV third generation sfermions and electroweak-inos [pj that have been proposed to 
ameliorate the SUSY flavour and CP problems are not in conflict with these LHC7 data. 

The LHC has performed spectacularly and has already delivered an integrated lumi- 
nosity of 5 fb -1 and, as we mentioned, is expected to deliver a comparable or larger data 
set in 2012. This motivated us to extend our earlier projections Q of the LHC7 reach for 
SUSY to integrated luminosities up to 30 fb -1 . As before, we work within the mSUGRA 
framework, the parameter space of which is given by, 

m , my 2 , A , tan /3, sign(p) . (1.1) 

Here, mo is a common GUT scale soft SUSY breaking (SSB) scalar mass, m 1 / 2 is a common 
GUT scale soft supersymmetry breaking (SSB) gaugino mass, Aq is a common GUT scale 
trilinear SSB term, tan j3 is the ratio of Higgs field vevs, and \i is the superpotential higgsino 
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mass term, whose magnitude, but not sign, is constrained by the electroweak symmetry 
breaking minimization conditions. 

For each model parameter space point, many simulated collider events are generated 
and compared against SM backgrounds with the same experimental signature ||]. A 6- 
dimensional grid of cuts[|| is then employed to enhance the SUSY signal over SM back- 
grounds, and the signal is deemed observable if it satisfies pre-selected criteria for observ- 
ability. Based on previous studies we include in our analysis the following channels: 

• jets + E™ 1SS (no isolated leptons), 

• U + jets + Ef iss , 

• two opposite-sign isolated leptons (OS) + jets + EJp lss , 

• two same-sign isolated leptons (SS)+jets + E™ 1SS , 

• U + jets + E^ iss . 

For the simulation of the background events, we use AlpGen|| to compute the hard 
scattering events and Pythia [n]] for the subsequent showering and hadronization. For the 
final states containing multiple jets (namely Z(—> 11, vv) + jets, W{— > lis) + jets, bb + jets, 
ti+jets, Z + bb+jets, Z+tt+jets, W+bb+jets, W+tt+jets and QCD), we use the MLM 
matching algorithm to avoid double counting. All the processes included in our analysis 
are shown in Table 1 of Ref. Q as well as their total cross-sections, number of events 
generated and event generator used. Here, we show in Table |3] the various backgrounds 
along with fc-factors 1 used to normalize the generator cross sections to NLO QCD results 
where available. The background fc-factors were computed using MCFM[jll|] for the NLO 
cross sections, and AlpGen for the LO ones. 

The signal events were generated using Isajet 7.79 (l^] which, given a mSUGRA pa- 
rameter set, generates all 2 — > 2 SUSY processes in the right proportion, and decays the 
sparticles to lighter sparticles using the appropriate branching ratios and decay matrix 
elements, until the parent sparticle cascade decay [p^j terminates in the stable lightest su- 
persymmetric particle (LSP), assumed here to be the lightest neutralino. Total gluino 
and squark production cross sections have been presented in Ref. J§ at NLO QCD using 



Prospino[14], and will not be repeated here. It is worth noting that for rriq ~ mg, gq 
associated production is the dominant strongly interacting SUSY production mechanism, 
while for rriq S> mg, gg pair production tends to dominate. 

For event generation, we use a toy detector simulation with calorimeter cell size Atj x 
Acp = 0.05 x 0.05 and —5 < ij < 5 . The HCAL (hadronic calorimetry) energy resolution is 
taken to be 80%/\/£e3% for \rj\ < 2.6 and FCAL (forward calorimetry) is W0%/VE®5% 
for |r/| > 2.6, where © denotes a combination in quadrature. The ECAL (electromagnetic 
calorimetry) energy resolution is assumed to be 3%/%/^© 0.5%. We use the cone-type 



x By fc-factor, here we actually mean a NLO /a LO . Normally, one compares the two cross sections using an 
identical renormalization/factorization scale for the two cases. Here, we merely compute a LO using AlpGen 
and a NLO using MCFM, using the pre-programmed default scale choices for the latter. 
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SM process 
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SM process 


/c-factor 


ti 


0.99 


QCD, bb 




Z/l + J ets 


1.47 


Z + ti 




W + jets 


1.53 


TJ7" I -4-T 
W + tt 




Z(-> vv) + bb 


1.18 


W + bb 




Z/rt-t 11) + bb 


1.03 


W + tb 




WW 


1.38 


titi 




wz 


1.47 


tibb 




zz 


1.35 


bbbb 





Table 1: Background processes included in this LHC7 study, along with the fc-factor (from MCFM 
and AlpGen) used (when available) to normalize to NLO QCD. For ti production, the renormaliza- 
tion scale is chosen to match the NLO cross section. The event generator used, total cross sections 
and number of generated events are listed in Table 1 of Rcf . || . All light (and b) partons in the 
final state are required to have Et > 40 GeV. For QCD, we generate the hardest final parton jet 
in distinct bins to get a better statistical representation of hard events. 



Isajet |Tj] jet-finding algorithm to group the hadronic final states into jets. Jets and 
isolated lepton are denned as follows: 

• Jets are hadronic clusters with \t]\ < 3.0, R = \J Arj 2 + Acfi 2 < 0.4 and Exijet) > 50 
GeV. 

• Electrons and muons are considered isolated if they have < 2.0, pr(l) > 10 GeV 
with visible activity within a cone of AR < 0.2 about the lepton direction, Y,E^ Us < 5 
GeV. 

• We identify hadronic clusters as 6-jets if they contain a B hadron with Et{B) > 15 
GeV, r]{B) < 3 and AR(B,jet) < 0.5. We assume a tagging efficiency of 60% and 
light quark and gluon jets can be mis-tagged as a 6-jet with a probability 1/150 for 
E T < 100 GeV, 1/50 for E T > 250 GeV, with a linear interpolation for 100 GeV 
<E T < 250 GeV |§. 

As in Ref. B, we define the signal to be observable if 



S > max 



5\fB, 5, 0.25 



where S and B are the expected number of signal and background events, respectively, 
for an assumed value of integrated luminosity. The requirement S > 0.2i? is imposed to 
avoid the possibility that a small signal on top of a large background could otherwise be 
regarded as statistically significant, but whose viability would require the background level 
to be known with exquisite precision in order to establish a discovery. For cases with very 
low signal and BG event numbers, we require the Poisson probability to correspond to the 
5cr level. 

The grid of cuts used in our optimized analysis is: 
• £?p iss > 50,100 - 1000 GeV (in steps of 100 GeV), 
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• n(jets) > 2, 3, 4, 5 or 6, 

• n(b — jets) > 0, 1, 2 or 3, 

• Et(Ji) > 50 - 300 GeV (in steps of 50 GeV) and 400-1000 GeV (in steps of 100 GeV) 
(jets are ordered j\ — j n , from highest to lowest Et), 

• E T (j 2 ) > 50 - 200 GeV (in steps of 30 GeV) and 300, 400, 500 GeV, 

• n(£) = 0, 1, 2, 3, OS, SS and inclusive channel: n(£) > 0. (Here, £ = e, //). 

• 10 GeV< m{£+£-) < 75 GeV or m{£+£-) > 105 GeV (for the OS, same flavor (SF) 
dileptons only), 

• transverse sphericity St > 0.2. 




500 1000 1500 2000 2500 3000 

m (GeV) 

Figure 1: The optimized SUSY reach of LHC7 for different integrated luminosities combining the 
different channels described in the text. The signal is observable if it falls below the solid contour for 
the corresponding integrated luminosity. The fixed mSUGRA parameters are Ao = 0, tan j3 — 45 
and /i > 0. Gluino mass contours are shown by the dashed, dark grey curves. The shaded grey area 
is excluded due to stau LSPs (left side of figure) or no electroweak symmetry breaking (right side 
of figure), while the shaded grey area marked "LEP excluded" is excluded by non-observation of a 
sparticle signal from LEP2 searches. All sparticle and background cross sections are normalized to 
NLO QCD values via /c-factors. 

We show in Fig. [l| the optimized 5a discovery reach of LHC7 for various choices of 
integrated luminosity in the mo vs. mpj plane. We also take Aq = 0, tan/3 = 45 and 
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\i > 0, with nit = 172.6 GeV. 2 Gluino iso-mass contours are shown, as obtained using the 
ISASUGRA routines[16] in Isajet. We see from Fig. p] that with ~ 1 ft -1 of integrated 
luminosity, the LHC7 sensitivity does indeed extend to m g ~ 1.1 TeV for niq ~ m g , and 
to nig ~ 0.65 TeV for nig 3> m g s . For 5 fb -1 of integrated luminosity (for which we expect 
ATLAS and CMS analyses in Spring 2012), the LHC discovery reach extends to m g ~ 1.3 
TeV for nig ~ nig, and to m g ~ 0.8 TeV for nig 3> m g . As integrated luminosity moves 
into the 20-30 ft -1 regime, the LHC7 reach for niq ~ m g moves up to m g ~ 1.5 — 1.6 
TeV. For the case where niq 3> m g , the 20-30 ft -1 LHC reach approaches m g ~ 1 TeV. 
We stress that - as discussed above - while non-observation of the signal at LHC7 may 
qualitatively point toward very heavy gluinos and first generation squarks, this does not in 
and of itself preclude SUSY as the new physics that stabilizes the weak scale [Q] because 
third generation squarks and electroweak-inos could still be at the sub- TeV scale. 

While our results are presented for the particular choice of mSUGRA parameters 
Aq = and tan f3 = 45, we emphasize here that we expect these results to hold also for 
other choices of Aq and tan/3, and also for n, < 0. Variation of Aq mainly affects third 
generation sparticle masses, while the reach is determined mostly by m g and the first 
generation squark masses. Moreover, variation of tan/3 mainly affects the size of b and r 
Yukawa couplings, and these feed only weakly into the reach plots: for instance, sparticle 
decays to third generation matter are enhanced at large tan/3[|2(| where 6-tagging may 



somewhat enhance the LHC reach for gluinos [21] as already demonstrated by ATLAS 

To give the reader an idea of the dominant event topologies in which experiments at 
LHC7 will be able to probe SUSY in the 2012 run, we show in Fig. [2] the optimized 5a 
reach via the 0£, 11, OS dilepton, SS dilepton and the trilepton channels for 20 fb . The 
striking feature of the figure is that while the reach is dominated by the low multiplicity 
(ni = 0, 1) lepton channels for niQ ~ 1.5 TeV, the reach in the low background but rate- 
limited trilepton channel becomes competitive with that in other channels if squarks are 
essentally decoupled at LHC7 as could well be the case. We have checked that this is true 
also for an integrated luminosity of 10 fb . 

In summary, we have presented updated 5a discovery contours for the paradigm 
mSUGRA/CMSSM SUSY model for LHC7 with 5-30 ft.- 1 of integrated luminosity. These 
results help to understand the capabilities of LHC7 for discovering supersymmetry in 2012- 
2013. Within mSUGRA, for integrated luminosity 20-30 fb -1 , we expect LHC7 to probe 
m g up to ~ 1.6 TeV for nig ~ m g , while we expect LHC7 to probe up to m g ~ 1 TeV for 
nig S> nig. If squarks are much heavier than gluinos, the reach at LHC7 via the inclusive 
trilepton channel will be competitive in reach with the canonical jets plus I?™ 1SS channel. 

2 Recent evidence from Atlas^j] and CMS[[^| using 5 fb -1 of data show some evidence for a Higgs scalar 
h with rrih ~ 125 GeV. For Ao = 0, it is very difficult to accommodate such a Higgs mass in the mSUGRA 
model. For Aq ~ ±2mo, then mj, ~ 125 GeV can be accommodated, but mainly at rather high mo ~ 2— 10 
TeV. For more details, see e.g. Ref. ps| . Our reach projections are largely insensitive to variation in Ao 
(and subsequent small changes in mt) as explained below. 

3 We stress that the curves presented here include an optimization over several search channels and 
correspond to a 5a discovery reach. Care must be taken when comparing these results with experimental 
bounds, which are usually presented for single channels at 95% C.L. (~ 2a). 
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Figure 2: The optimized ha SUSY reach of LHC7 in various channels classified by lepton multi- 
plicity: 0£, 1£, SS dilepton, OS dilepton and trilepton for an integrated luminosity of 20 fb _1 . Any 
mSUGRA point will be observable if it falls below the corresponding contour. The fixed mSUGRA 
parameters are Aq = 0, tan/3 = 45 and /i > 0. Gluino mass contours are shown by the dashed, 
dark grey curves. The shaded grey area is excluded due to stau LSPs (left side of figure) or no 
electroweak symmetry breaking (right side of figure), while the shaded grey area marked "LEP 
excluded" is excluded by non-observation of a sparticle signal from LEP2 searches. All sparticlc 
and background cross sections are normalized to NLO QCD values via fc-factors. 
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